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We investigate the low-energy properties of the orthogonal-dimer spin chain characterized by a 
frustrated dimer-plaquette structure. When the competing antiferromagnetic couplings are varied, 
the first-order quantum phase transition occurs between the dimer and the plaquette phases, which 
is accompanied by nontrivial features due to frustration: besides the discontinuity in the lowest 
excitation gap at the transition point, a sharp level-crossing occurs for the spectrum in the plaquette 
phase. We further reveal that the plateau in the magnetization curve at 1/4 of the full moment 
dramatically changes its character in the vicinity of the critical point. It is argued that the first- 
order phase transition in this system captures some essential properties found in the two-dimensional 
orthogonal-dimer model proposed for SrCu2(BOa)2. 

PACS numbers: 75.10Jm, 75.40Cx 
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I. INTRODUCTION 

Recent extensive experimental and theoretical inves- 
tigations on low-dimensional spin systems with frustra- 
tion have been providing a variety of interesting topics. 
Among others, the -spin-gap compound found recently 
by Kageyama et alE, SrCu2 (B03)2, exhibits a number 
of nontrivial properties due to strong frustration.ETLj In 
particular, the discovery of the magnetization plateaus 
at 1/3, 1/4 and 1/8 of the full momcaitUtl has been 
stimulating further intensive studiesjjcl jlhe remark- 
able point claimed by Miyahara and Uedacl is that this 
frustrated system is a prototypical example of the two* 
dimensional (2D) version of the orthogonal-dimer models 
whose uniajjc— structure gives rise to various unusual 
properties£3i!3 This 2D model is known to be equiv- 
alent to the frustrated Shastry-Sutherland. model on a 
square lattice with some diagonal bonds O It has been 
recently shown that there exists a novel first-order quan- 
tum phase transition between the dimer and plaquette 
phases in the 2D orthogonal-dimer model, which is ac- 
companied by the jump and Cjusp singularities in the spin 
gap near the transition point .113 What is most interesting 
is that the compound SrCu 2 (603)2 may beJacated in 
the vicinity of the first-order transition point nu 

As seen from the above studies, characteristic proper- 
ties in the 2D system around the first-order transition 
point are certainly caused by the strong frustration com- 
mon to this class of the orthogonal-dimer spin systems. 
Therefore, to clarify the essential properties of the sys- 
tem, further systematic studies on the quantum phase 
transition are highly desired. The ID version of the 
orthogonal-dimer modecjlij may be the most appropri- 
ate system for this purpose, because it is the simplest 
model which possesses the frustration effect due to the 
dimer-plaquette structure. Motivated by these hot top- 
ics, in this paper, we study the first-order transition in 
the ID orthogonal-dimer spin chain in detail£3~ll3 The 
Hamiltonian we shall deal with is. 



FIG. 1. The orthogonal-dimer chain. The solid and dashed 
lines indicate the antiferromagnetic exchange couplings J and 
J', respectively. 
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where S, is the s = 1/2 spin operator at the i-th site, and 
the indices and < i,j > represent the summation 

over intra- and inter-dimer pairs, respectively (see Fig. 
[j]). The magnetic field is denoted as H for which we set 
5/iB = 1 for convenience. Both of the exchange couplings 
J and J' are assumed to be antiferromagnetic. We shall 
use the normalized parameters j = J' / J and h = H/J 
in the following discussions. .— . 

It was previously shown by Ivanov and RichterliZI that 
the above orthogonal-dimer chain undergoes the first- 
order quantum phase transition between the dimer phase 
(for small j) and the plaquette phase (for large j) as the 
coupling ratio j is varied. However, the unique proper- 
ties inherent in this frustrated system have not been dis- 
cussed in detail yet, especially around the critical point. 
In what follows, by means of the exact diagonalization 
(ED) and series expansion methods, we demonstrate that 
the excitation gap in this system exhibits nontrivial be- 
haviors such as the jump and cusp singularities around 
the transition point, reflecting the dimer-plaquette dual 
properties characteristic of this frustrated system. We 
also point out that such properties are not specific to 
the ID system but also common to the 2D orthogonal- 
dimer system. By computing the magnetization curve by 
means of the ED together with the density-matrix renor- 
malization group (DMRG),li3E3 we further reveal that 
the formation of the magnetic plateaus is dramatically 
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affected by the dimer-plaquette dual properties, and the 
resulting magnetic phase diagram has a rich structure. 

This paper is organized as follows. In Sec. |l[ we in- 
vestigate how the characteristic properties inherent in the 
orthogonal-dimer chain emerge in the low-energy excita- 
tions around the first-o rcLer .transition point by means 
of the ED, the DMRGB'B and the series expansion 



method In Sec. |IE, by calculating the magnetiza- 



tion curve by the ED and the DMRG, we show that the 
plateau formation in the magnetization changes its char- 
acter around the transition point, reflecting the strong 
frustration effects. The last section is devoted to sum- 
mary and discussions. 



II. ZERO-FIELD PROPERTIES 

In this section, we investigate the quantum phase tran- 
sition in the orihogonal-dimer chain in the absence of the 
magnetic fieldJiJ bv exploiting the ED and the series ex- 
pansion methods E3 



A. Ground state 

We start with the ground state properties in the dimer 
phase. When j — 0, the system is reduced to an as- 
sembly of the decoupled dimers denoted by the solid line 
in Fig. [I], for which the product of independent dimer- 
singlets gives the ground state. The remarkable point for 
the orthogonal-dimer chainBO is that this simple dimcr- 
singlet state is always an exact eigenstate of the Hamil- 
tonian ([!]) with the energy E g /JN = —3/8 in the entire 
range of j (see Fig. |[), where N is the number of total 
sites. Accordingly, the dimer-singlet state should be the 
exact ground state up to a certain critical value of j. 
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FIG. 2. The ground state energy per site for the orthog- 
onal-dimer chain as a function of j{= J'/J)- The flat line 
(E S /JN = —3/8) indicates the exact ground state energy of 
the dimer phase. The solid line is the result obtained by the 
plaquette expansion. We also show the results for the finite 
chain of N = 24 as open circles. 



In the opposite limit of large j, the ground state is 
given by the disordered singlet state which is adiabati- 
cally connected to the isolated plaquette-singlets denoted 
by the broken line in Fig. [j]. Therefore, starting from 
the isolated plaquette singlets for j = oo, we can evalu- 
ate the ground state energy E g for the plaquette phase 
with finite j by means of the series expansion. Perform- 
ing the plaquette expansion up to the eleventh order in 
j^ 1 combined jSidth the first-order inhomogeneous differ- 
ential methodj£3 we have obtained the ground state en- 
ergy E g rather precisely for the plaquette phase. The 
result is shown as the solid line in Fig. ||. It is seen that 
the ground state energy of the plaquette state coincides 
with that of the dimer state at a certain value of j, at 
which the first-order quantum phase transition occurs. 
The critical value is estimated as j c — 0.8f900, which 
is further confirmed to be accurate up to the above fig- 
ure by the DMRG calculation for the infinite chain. We 
note here that the present results for the ground state ate 
consistent with those obtained by Ivanov and Richter.liZI 



B. Spin excitations 

Let us move to the spin excitation spectrum. In the 
dimer phase (0 < j < j c ), we can construct a low-energy 
triplet excitation by substituting a local triplet for one of 
singlet-dimers forming the ground state. The remarkable 
point is that the hopping of this local triplet across the 
singlet-dimer such as the f-3 dimer in-Fig. I is forbidden 
by its characteristic crystal structure.u This implies that 
an excited triplet state is completely localized inside the 
finite strip (N = 6), and has no dispersion for its spec- 
trum. Thus we can exactly estimate the spin gap from 
the ED calculation for the cluster of N = 6. In Fig. ||, 
the ED results are shown as the solid line in the dimer 
phase. 
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FIG. 3. The spin gap A/J for the orthogonal-dimer chain 
as a function of j(= J'/J)- The crosses and open circles 
indicate the results obtained by the ED method for the finite 
chain of N = f 6 and 24, respectively. See the text about the 
solid and broken lines. 

In the plaquette phase (j > j c ), it is naively expected 
that the low-energy excitation is described in terms of a 
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triplet excitation on the isolated plaquette. However, the 
frustrating diagonal interaction in each plaquette should 
considerably affect the excitation spectrum. To see this, 
it is instructive to examine the energy-level structure of 
the isolated plaquette with the diagonal bond (see Fig. 
|l|), whose Hamiltonian reads 

Hplaquottc = JSl • S3 + J' (si + S3) • (s 2 + S 4 ) . (2) 

We list the energy eigenvalues of the plaquette in Table Q, 
where S13 (S24) is a quantum number of the spin Si + S3 
(S2 + S4) and S is that for the total spin. It is seen in 



TABLE I. Eigenstates of an isolated plaquette. 
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this table that the eigenstates of the isolated plaquette 
are classified into two sectors: the D-sector for S13 = 
and the P-sector for Si 3 = 1. For j > 1/2, the ground 
state is always singlet in the P-sector, whereas both of 
two separated sectors appear for the excitations when the 
coupling ratio j is varied. When 1/2 < j < 1, the lowest 
excitation is given by the S24-singlet and S24-triplet in 
the D-sector, which are four-fold degenerate with the ex- 
citation energy AE/J — 2j — 1. In the regime j > 1, the 
first excitation is a P-sector triplet with the excitation 
energy AE/J = j. Note that as far as the singlet-triplet 
excitations in the P-sector are concerned, there are no 
contributions from the D-sector. This simple fact allows 
us to replace S13 in the plaquette with the s = 1 spin, 
which will be used below to perform the mixed-spin clus- 
ter expansion for the P-sector. 

Keeping the above properties in mind, we now turn 
back to the orthogonal-dimer chain in the thermody- 
namic limit, by taking into account the interaction among 
plaquettes. Note that excitations in the D-sector are 
still completely separated from those in the P-sector, ac- 
cording to the constraint specific to the orthogonal-dimer 
structure. This remarkable fact enables us to study the 
excitation spectrum for the D-sector and the P-sector 
separately. Namely, we perform the plaquette expan- 
sion up to the eleventh order in j -1 for the D-sector, 
whereas for the, P-sector we exploit the mixed-spin clus- 
ter expansions up to the seventh ordepjn j _1 with a 
starting spin configuration 1/2 o 1 o 1/2. Applying the 
first-order inhomogeneous differential methodc3 to the 
obtained series, we then deduce the lowest excitations 
both for the D-sector and the P-sector, as shown in Fig. 
|3|. In this figure, the dashed line in the plaquette phase 
represents the energy for four-fold degenerate excitations 
in the D-sector whereas the solid line is for a triplet ex- 
citation in the P-sector. It is seen that the energy of two 
kinds of excitations intersects each other at j' c — 0.872, 



giving rise to a cusp singularity in the spin gap as a func- 
tion of j. In the region j > j' c , both of the ground state 
and the lowest excited state belong to the P-sector on 
each plaquette, as is the case for the simple plaquette 
chain.lij On the other hand, when j c < j < j' c , the low- 
est excitation is described by four-fold degenerate level 
in the D-sector although the ground state still belongs 
to the plaquette phase (P-sector). The present plaque- 
tte expansion further uncovers that the wave function 
of the localized D-sector excitation with no dispersion is 
spatially extended over a number of sites, which shows 
sharp contrast to the dimer phase where the wave func- 
tion of the triplet excitation is completely localized at a 
given site. This may imply that the four-fold degenerate 
excitations characterized by the D-sector possess the in- 
termediate properties between those typical for the dimer 
phase and the plaquette phase. As j is further decreased, 
we encounter the first-order quantum phase transition, at 
which the ground-energy of the two phases coincide with 
each other, while the first excited states have still differ- 
ent energies. Therefore, the spin gap at the transition 
point jumps from A + (j c ) = 0.3590(2) for the plaquette 
phase to A~(j c ) = 0.32309 for the dimer phase. 

In order to complement the series-expansion results, 
we also compute the lowest triplet excitation by the ED 
method for the finite chain (N — 16, 24) with the periodic 
boundary condition. The obtained spin gap A is shown 
as the crosses and the open circles in Fig. [|. It is seen 
that they are in fairly good agreement with the results 
obtained by the series expansion method, except for the 
vicinity of the first-order transition point j c , where the 
finite-size effect in the ED still remains. Here it should 
be noticed that the present results provide much more 
detailed information for spin excitations than those ob- 
tained by Ivanov and Richter, where the spin gap as a 
function o£-d has only a cusp-like singularity at the crit- 
ical point .E_ZJ We also note that the ladder system with a 
similar orthogonal-dimer structure-has the jump and the 
cusp singularities in the spin gapEJ 

III. PLATEAUS IN THE MAGNETIZATION 
CURVE 

In this section, we study the magnetic properties of 
the orthogonal-dimer chain. By using the ED and the 
DMRGEiiEj methods, we calculate the magnetization 
curve as a function of the magnetic field h(— H / J), and 
show that their characteristic behavior stems from the 
dimer-plaquette dual properties found for the excitations 
in the previous section. To clearly understand the magne- 
tization curve in the low-field region, we shall use the de- 
scription based on hard-core bosons, in which low-energy 
triplet excitations are regarded as hard-core bosons. 
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A. Dimer phase 

Let us start with the dimer phase. In Fig. ||, we show 
the magnetization curve for j — 0.7(< j c ) obtained by 
the ED method for the finite chain of N = 8, 16 and 24. 
We also show the results obtained by the DMRG method. 
In this figure, the plateaus appear in the magnetization 
curve at 1/4 and 1/2 of the full moment (to = 1/2) 
clearly. The mechanism of the 1/2-plateau formation is 




FIG. 4. The magnetization m(= M/N) as a function of h 
in the dimer phase (j = 0.7). Three thin lines represent the 
magnetization curves for the finite chain of N = 8, 16 and 24. 
The bold line is the result obtained by the DMRG method 
for the infinite chain. 

simply understood according to the unit-cell structure of 
the chain. As will be shown below, the 1/2-plateau ap- 
pears in the entire range of j. On the other hand, the 
1/4-plateau may be realized by the commensurate state 
in which each triplet excitation is aligned periodically at 
every finite- strip of N = 8, reflecting the completely- 
localized nature of the triplet excitation discussed in the 
previous section. Since the localized triplet does not cor- 
relate with those in the neighbor strips, the macroscopic 
number of the triplets become degenerate at the lower- 
critical field h c (= the spin gap A), where the magnetic 
first-order transition occurs. In Fig. |], we can indeed 
see that the magnetization curve jumps from m — to 
rn = 1/8 (1/4-plateau) at h c without any finite-size cor- 
rection. Increasing the magnetic field, the other first- 
order transition to the 1/2-plateau state occurs. It is 
seen in Fig. [| that there still remains the large-finite size 
effect near this critical point, since the triplet excitations 
correlate with each other in contrast to the transition 
discussed above. 

We here note that several different phases observed in 
the magnetization process are specified by the spatial ar- 
rangement of the spin quantum number, S13, for the di- 
agonal bond in each plaquette (see Table |). In the dimer 
phase, the P-sector (Si3=0) is realized in each plaquette 
and in the 1/4-plateau state, the D- and the P-sectors are 
crystallized alternately. On the other hand, in the region 
of m > 1/4, each diagonal bond in the plaquette forms 
S13 = 1 (the P-sector), and hence the frustration does 



not play any role for the magnetization process in this re- 
gion. Therefore by calculating the magnetization for the 
unfrustrated plaquette chain, we can simply reproduce 
that for the orthogonal-dimer chain in this region. 



B. Plaquette phase 

We recall that there are two sort of excitations char- 
acterized by the P- and the P-sector in the plaquette 
phase, which may provide a variety of the magnetization 
process. More precisely, the hard-core boson description 
suggests three possible behaviors for the magnetization 
curves: 

• case I: P-triplet excitations lie energetically lower 
than P-triplct excitations. The magnetization 
curve has a structure similar to that in the dimer- 
phase. 

• case II: P-triplet excitations lie slightly above the 
bottom of the P-triplct dispersion curve (see Fig. 
^j. Thus both of the features specific to the P- 
triplct and the P-triplet excitations appear in the 
magnetization curve. 

• case III: P-triplet excitations lie sufficiently higher 
above the spin gap for P-triplet excitations. The 
magnetization curve may be characterized solely by 
the P-triplet excitation. 

Since the above classification is simply based on the hard- 
core boson description, the interaction among bosons be- 
comes more important when the number of bosons (i.e. 
the strength of the external field) is increased, which re- 
quires more proper discussions beyond the hard-core bo- 
son description. In what follows, based on the numerical 
results, we confirm that the magnetization curves are in- 
deed classified into the above three categories. 

We start with the analysis of the case I. We show the 
magnetization process for j = 0.86 in Fig. [5], where we 
can see the 1/4- and 1/2-plateaus clearly. In particular 
the magnetization jumps from zero to the 1/4 plateau 
at the critical field h c with a small finite-size correction. 
In the plaquette phase, the 1/4-plateau is generated by 
the commensurate state in which the plaquette singlet 
and triplet in the P-sector are crystallized alternately. 
The characteristic point in the 1/4-plateau state is that 
P-triplet excitations interact with each other when they 
are condensed to form the 1/4-plateau, in contrast to 
those in the dimer phase where each triplet excitation is 
localized completely. We can indeed see this interaction 
effect in the ED spectrum, where the triplets gain the 
condensation energy, by which the critical field h c shifts 
to a slightly lower field than the spin-gap value. We note 
here that since the mechanism of the 1 / 4-plateau forma- 
tion (first-order transition) at h c is essentially the same 
in the dimer phase and in the plaquette phase, so that 
the critical field h c does not show the discontinuity at 
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FIG. 5. The magnetization m{= M/N) as a function of h 
in the plaquette phase (j = 0.86). Three thin lines represent 
the magnetization curves for the finite chain of N = 8, 16 and 
24. The bold line is the result obtained by the DMRG for the 
infinite chain. 



tinuously and jump to the 1/4-plateau. Beyond the 1/4- 
plateau, we again encounter the first-order phase transi- 
tion, which is accompanied by the jump in m. This jump 
singularity is followed by the continuous increase in to, 
as already discussed in the case I. 
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j = j c as a function of j (see Fig.gJ) , being contrasted 
to the discontinuity observed in the spin gap at j = j c . 
Increasing the field, the magnetic first-order transition 
takes place to the state in which S13 = 1 is formed in 
each plaquette. Beyond this critical field, the magne- 
tization traces the curve which is the same as that for 
unfrustrated plaquette chain. 

In the case II, the D-triplet excitation level lies slightly 
above the bottom of the P-triplet dispersion curve, as 
shown in Fig.^. The coexistence of two-kinds of distinct 




FIG. 6. The dispersion relation for j = 0.94. The flat 
line and the solid curve represent the four-fold degenerate 
D-triplet excitations and the P-triplet excitations, which are 
obtained by the mixed-spin cluster expansion and the plaque- 
tte expansion, respectively. 

excitations influences the magnetization curve consider- 
ably, as seen in Fig. ^. As h increases beyond the critical 
field (h c = A), the magnetization should develop with 
m ~ [h — he) 1 ' 2 , since it is dominated by the P-triplct 
excitation whose dispersion relation is quadratic near the 
bottom (q = tt/2). In fact, we can see a staircase struc- 
ture below the 1 /4-plateau for the ED results in Fig. |t], 
implying that the jump is now changed to the continuous 
increase in the magnetization. With slightly increasing 
h, however, the magnetization may stop to increase con- 



FIG. 7. The magnetization m(= M/N) as a function of h 
in the plaquette phase (j = 0.94). Three thin lines represent 
the magnetization curves for the finite chain of N = 8, 16 and 
24. The bold line is the result obtained by the DMRG for the 
infinite chain. 

In the case III, the excitation energy of the -D-triplet is 
pushed up in higher energy region of the dispersion curve. 
Thus, before the Zeeman energy lowers the D-triplet level 
down to zero, a number of bosons are accommodated 
in the P-triplct excited levels, and then the interaction 
effect between bosons becomes too significant to use the 
hard-core boson description. In this parameter region, 
we find that the 1/4-plateau completely disappears and 
the magnetization curve becomes smooth up to the 1/2- 
plateau. 




FIG. 8. The magnetization m(= M/N) as a function of h 
in the plaquette phase (j = 1.1). Three thin lines represent 
the magnetization curves for the finite chain TV = 8, 16 and 
24. The bold line is the result obtained by the DMRG for 
the infinite chain. Note that the 1/2-plateau state still exists 
near h = 1.9. 

Consequently, we end up with the magnetic phase dia- 
gram on the (j — h) plane as shown in Fig. ^. In this fig- 
ure, (a) and (a') indicate the dimer-singlet and plaquctte- 
singlet phase, respectively. The region (b) exhibits the 
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FIG. 9. The ground-state phase diagram on the j versus 
h plane. The phase boundaries are determined by the ED 
method for the finite chain of TV = 24 and the DMRG method 
for the infinite chain. 

plateau in the magnetization curve at 1/4 of the full mo- 
ment. We can see that the 1/4-platcau phase with the 
dimer-singlet structure exists not only in the dimer phase 
(j < jc) but also in a part of the plaquette phase. In the 
region (c), the 1/2-plateau appears, which is stable in the 
whole parameter region of j. All spins are polarized in 
(d). The bold (solid) line shows the phase boundary of 
the first- (second-) order transition. It is now seen that 
the magnetic phase diagram possesses quite a rich struc- 
ture, which is indeed caused by the dimer-plaquette dual 
properties inherent in this strongly frustrated spin chain. 
Here we wish to mention that the .ladder system with a 
similar orthogonal-dimer structural has the jurtm and 
the cusp singularities in the magnetization curve.t3 

IV. SUMMARY AND DISCUSSIONS 



of the full moment, according to how two kinds of triplet 
excitations change their relative positions energetically. 
We have thus found that the magnetic phase diagram on 
the (j — h) plane has a rich structure. 

Although we have dealt with the one-dimensional chain 
in this paper, we believe that the present study cap- 
tures some essential properties of the first-order phase 
transition and the resulting rich structure in low-lying 
excited states even in two dimension. For the com- 
pound SrCu 2 (B0 3 ) 2 , it is indeed pointed out that the 
low-energy triplet excitation has quite small mobility re- 
flecting the orthogonal-dimer structure.^! Since this com- 
pound is located in the vicinity of the first-order phase 
transition point, the present results encourage us to in- 
terpret the plateau-formation mechanism of the 2D com- 
pound in terms of the dual properties originating from 
the dimer and plaquette structures. In particular, it is an 
interesting problem to study how the characteristic prop- 
erties in ID are changed into those in 2D by adiabatically 
introducing the inter-chain coupling, which may provide 
a clue to fully understand the physics in 2D orthogonal- 
dimer system. 
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We have studied the low-energy properties of the 
orthogonal-dimer spin chain with a frustrated dimer- 
plaquette structure, by means of the exact diagonaliza- 
tion, the density matrix renormalization group and the 
series expansion method. When the inter-dimer cou- 
plings are varied, the first-order quantum phase tran- 
sition occurs at the critical value j c = 0.81900 between 
the dimer-singlet and plaquette-singlet phases. In the 
dimer phase, the ground state is exactly given by the de- 
coupled dimer-singlet state, and triplet excitations over 
it are completely localized. In the plaquette phase, on 
the other hand, we have found that the excited states 
are characterized by the dual properties: the coexistence 
of plaquette- like and dimer- like excitations. As a result, 
both of the cusp and the jump appear in the spin gap as 
a function of j, which are caused by the level-crossing be- 
tween the excited states. We have also clarified how the 
above dual properties influence the magnetization pro- 
cess of the chain by analyzing the numerical results in 
terms of the hard-core boson description. In particu- 
lar, we have shown that three types of the magnetization 
curves appear around the magnetization plateau at 1/4 
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